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We discuss the chiral phase diagram in the parameter space of lattice QCD with 
minimal-doubling fermions, which can be seen as lattice fermions with flavored 
chemical potential terms. We study strong-coupling lattice QCD with the Karsten- 
Wilczek formulation, which has one relevant parameter ^ as well as gauge coupling 
and a mass parameter. We find a nontrivial chiral phase structure with a second- 
order phase transition between chiral symmetric and broken phases. To capture 
the whole structure of the phase diagram, we study the related lattice Gross-Neveu 
model. The result indicates that the chiral phase transition also exists in the weak- 
coupling region. From these results we speculate on the chiral phase diagram 
in lattice QCD with minimal-doubling fermions, and discuss their application to 
numerical simulations. 
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I. INTRODUCTION 



The famous no-go theorem [l, 2] states that lattice fermion actions with chiral symmetry, 
locality and other common features must produce massless degrees of freedom in multiples 
of two in a continuum limit. This is contradictory with a phenomenological fact that there 
exist only three quarks with masses below the QCD scale. By now several fermion construc- 
tions to bypass the no-go theorem have been developed, although all of them have their 
individual shortcomings: The explicit chiral symmetry breaking with the Wilson fermion 
approach results in an additive mass renormalization, which in turn requires a fine-tuning 
of the mass parameter for QCD simulations. Domain- wall 4,5] and overlap fermions 0, Q] 
produce a single fermion mode by modifying the definition of chiral symmetry or introducing 
the momentum-dependent chiral charge, although they lead to rather expensive simulations 
algorithms.. These approaches attempt to realize single fermionic degrees of freedom by 
breaking the requisite conditions for the no-go theorem. On the other hand, there is another 
direction to approach numerical simulations. According to l|, Hypercubic symmetry and 
reflection positivity of actions result in 2 d species of fermions where d stands for the dimen- 
sion. Thus it is potentially possible to reduce the number of species by breaking hypercubic 
symmetry properly. Actually, the staggered fermion approach 8KL2J, with only 4 species 
of fermions does this and possesses flavored-hypercubic symmetry instead. However this 
requires rooting procedures for the physical 2 or (2 + l)-flavor QCD simulation. 

A possible goal in this direction is a lattice fermion with 2 species, the minimal num- 
ber required by the no-go theorem. Such a minimal-doubling action was first proposed by 



Karsten, and later by Wilczek 
known as Creutz-Borici type 



la-hjz 



Other than the original t ype , two more types are 
17| and Twisted-ordering type 18j|. These fermions all 



possess one exact chiral symmetry and ultra locality. As such they could be faster for simu- 
lation, at least for two-flavor QCD, than other chirally symmetric lattice fermions. However 
it has been shown 19M22| that we need to fine-tune several parameters for a continuum limit 
with these actions. This is because they lack sufficient discrete symmetry to prohibit redun- 
dant operators from being generated through loop corrections 23l-l26|. Thus the minimally 
doubled fermions have not been extensively used so far. Nevertheless, there is the possibil- 
ity to apply them to simulations if one can efficiently perform the necessary fine-tuning of 
parameters. 
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In this paper we pursue the chiral phase structure in the parameter plane for minimal- 
doubling lattice QCD for two purposes: One purpose is to understand properties of these 
formulations theoretically. The other purpose is to show their applicability to lattice QCD 
simulations since understanding a nontrivial phase diagram, as with the Aoki phase in 
Wilson fermion 27H31|. can be useful in showing the applicability to lattice QCD. We first 
show that minimal-doubling fermions can be seen as a special case of lattice fermions with 
species-dependent (imaginary) chemical potential terms. We focus on the Karsten-Wilczek 
formulation with one relevant parameter /x 3 , which corresponds to a mass parameter in the 
analogy of Wilson fermion. We next investigate the chiral phase structure in the space of 
the gauge coupling and this relevant parameter. We analyze strong-coupling lattice QCD, 
and show that chiral symmetry is spontaneously broken in a certain range of the parameter 
while the chiral condensate is zero outside the range. We find a second-order phase transition 
between these chiral symmetric and broken phases. We also show that pion becomes massless 
as a Nambu-Goldstone boson in the chiral-broken phase while the sigma meson becomes 
massless on the second-order phase boundary due to the critical behavior of the second- 
order phase transition. We investigate the lattice Gross-Neveu model to capture an entire 
phase structure. From these results we suppose a similar chiral phase structure in 4d QCD 
with flavored-chemical-potential lattice fermions, and discuss their applicability to lattice 
QCD simulations. 

In Sec. Ull we study lattice fermions with flavored chemical potential, or minimal-doubling 
fermions. In Sec. 1111} we investigate a chiral phase structure in the framework of strong- 
coupling lattice QCD. In Sec. IIVI we study the Gross-Neveu model to obtain information 
of the whole phase diagram. In Sec. |Vj we discuss a phase structure in 4d QCD from the 
results of the last two sections. Section [VI] is devoted to a summary and discussion. 



II. FLAVORED-CHEMICAL-POTENTIAL FERMION 

In this section we study minimal-doubling fermions and their generalization. Before 
going to the main theme, let us remind ourselves of Wilson's prescription to shirk the no-go 
theorem. Wilson fermion extracts one light fermion by introducing a species-dependent mass 
term, which we call a "flavored-mass term" . A free action and a Dirac operator of Wilson 



i 



fermion are given by 



Sw = a ^ ~ + a }_^^n ^2 ' W 



aD w (p) = r^sinp^a + ^(1 - cos p^a). (2) 

We here exhibits a lattice spacing a to manifest mass dimensions of each term. In this 
formulation, 15 out of 16 species have 0(l/a) mass and are decoupled in the naive continuum 
limit. As shown in 32|, this is not the only case of flavored-mass terms. There are four types 
of nontrivial flavored-mass terms (Mf = Mp, My, Ma, Mr), which satisfy 75 hermiticity and 
give second-derivative terms up to 0(a 2 ) errors as with the usual Wilson term. A general 
form of Wilson-type fermions are written as 

aDfJjp) = ry M sinp M a + M { (p), (3) 

where Mf(p) stands for flavored-mass terms. Details of species-splitting depends on explicit 
forms of Mf(p) js^] . 

We now consider further deformation of the fermion action. We multiply Mf by 74 or ^74 

as 

aD ic (p) = r^sinp^a + (i)j 4 M f (p). (4) 

In these cases, degeneracy of species is lifted by specie-dependent real or imaginary chemical 
potential terms, not by species-dependent mass. It means that we can get rid of some 
doublers by this method too. We name such terms as "flavored-chemical-potential terms" , 
and name lattice fermions with them as "flavored-chemical-potential (FCP) fermions" . (We 



will later discuss problems with this kind of naive introduction of chemical potential |37|.) 
It is obvious that a real type of FCP terms breaks down 75 hermiticity and leads to a sign 
problem while an imaginary type keeps it and has no sign problem. An outstanding point 
in this formulation is that ultra- locality and one exact U(l) chiral symmetry remains intact 



as 



{ 75 ,Ac(p)} = 0. (5) 



(See Refs. [32j-[36j for details of chiral symmetry in this type of lattice fermions.) In principle, 
by using this deformation we can reduce 16 species to smaller multiple numbers of two 
without losing all chiral symmetries. We note that the chemical potential term, of course, 
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FIG. 1: Species-splitting in Wilson and Karsten-Wilczek fermion. Circled numbers stand for the 
number of massless flavors on each point. 
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breaks hypercubic symmetry into cubic symmetry and breaks C, P, T into CT and P 
as far as Mf(p) is cubic-symmetric. It means that this formulation automatically corresponds 
to finite-density systems unless we tune several parameters. In this paper we concentrate 
on the following explicit form of imaginary-type FCP fermions, 

3 



n,fi 



2a +a 



^2 ^«* 74 



2a 2 



(6) 



aD KW (p) = ij^sinp^a + ry 4 ^(l 

3=1 



COS PjCl). 



(?) 

Here 14 species is decoupled in the naive continuum limit while two species at p = (0, 0, 0, 0) 
and p = (0,0, 0,7r/a) has zero mass and zero imaginary chemical potential 52[. More 
precisely, among 16 species, two species have zero imaginary chemical potential, six have 
2/ a, six have 4/a and two have 6/a. In Fig. [T]we compares specie-splitting of KW fermion 
in chemical potential space to that of Wilson fermion in mass space. It is notable that 



two-flavor is the minimal number allowed 

ikfKW) fermion 



as the Karsten-Wilcze 



byt 
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le no-go theorem. This form has been known 



14J, which is the first known type of "minimal- 



doubling fermions" 13l4l8|. It has one exact chiral symmetry, ultra- locality, cubic symmetry, 
CT and P. Since the chemical potential term breaks discrete symmetries into the subgroup, 
we need to fine-tune three parameters for one dimension-3 (^74^) and two dimension-4 
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(ip^fzdiip, Fj^Fj^) counterterms in order to ta' 



the zero(T,/x) lattice QCD simulations 24|-|26|. 



ie a Lorentz-symmetric continuum limit for 



Among the three counterterms, in this paper we mainly deal with the dimension-3 term 
fistpi^ip with a relevant parameter /i 3 and the dimension-4 term d^^d^ip with a marginal 
parameter since we study the strong-coupling lattice QCD and the Gross-Neveu model, 
which contain no plaquette action. In particular the parameter /i 3 is of special importance: 
It changes the number of flavors and plays an important role in the chiral phase structure. 
Furthermore the quantum effects produce 0(1 /a) additive chemical potential renormaliza- 
tion in this case instead the additive mass renormalization, and we need to cancel it by 
adjusting // 3 even for the application to the imaginary-chemical-potential lattice QCD. This 
necessity of parameter tuning is also understood from the well-known fact that the naive 
introduction of chemical potential into lattice fermions leads to divergence of energy density 
and requires a counterterm due to the violation of the abelian gauge invariance as shown in 



Ref. [37]. 



We here write the KW fermion action of the interacting theory as 



<Skw — ^ 



1 4 



2 



3 

3=1 

+ y$s74 (Un,n+4*Pn+4 ~ £4,n-4^n-4) 



(8) 



where we introduce a parameter r in analogy with the Wilson parameter. We introduce the 

dimension-3 counterterm with the parameter /x 3 . Although we mainly focus on fi 3 in this 

work, we also introduce the dimension-4 counterterm with the parameter d^ in Eq. (JS} which 

is also relevant for the strong-coupling study. We make all the quantities dimensionless. Now 

let us look into how the number of flavors depend on /i 3 . For a free theory, the associated 

massless Dirac operator in momentum space is 

4 3 
clDkw(p) = i 7^ sin ap^ + z74(/i 3 + 3r — r cos apj + <i 4 sin ap±). (9) 
n=\ j=l 

We first look into a minimal-doubling range for /x 3 and the speed of light in the range. We 
for simplicity take r = 1. For — 1 — d^ < /i 3 < 1 + d±, we have only two zeros, both of which 
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have the form as p = (0, 0, 0, ^ arcsin(— j^-))- By expanding the momentum as p = p + q 
around the zeros, the coefficient of in the Dirac operator is given as 

3 

(1 + d 4 ) [sin aj?4 • cos ag 4 + cos ap^ ■ sin aq^] + /i3 + 3 — >^ cos(p; + qi) 



i=i 



= [(1 + gLl) cosap4] a^4 + 0(a 2 q 2 ). (10) 



Therefore, for general values of us in this range, the speed of light is modified as 



g 4 (l + d 4 )\ 1 



(1 + ^ 4 ) 2 




+ 0{aq 2 ). 



(11) 



We can fix this speed of light to a correct value up to the 0(a) discretization errors by taking 
(1 + d^) 2 = 1 + We will discuss details of Lorentz symmetry restoration in Sec]Vj Now, 
let us classify fi% parameter regions by the number of physical flavors in the free theory. 
Since d± just gives a shift of the parameter regions, we now consider d± = for simplicity. 
The above minimal-doubling range is given just by — 1 < ^3 < 1. For — 7 < 113 < — 5, we 
again have two zeros. For fi 3 < —7 and fi^ > 1, there is no zero of the Dirac operator, which 
means that there are no physical fermions. For —5 < ^3 < —3 and —3 < fi^ < —1, six 
zeros exist. [13 = 1, —7 have one zero, but the dispersion relations have a unphysical form 
as ~ p + p\. [13 = — 1, —5 have four zeros, but the dispersions are unphysical. fi^ = —3 has 
six zeros, but the dispersion is again unphysical. In the end, the parameter ranges where 
physical fermions can be described are —7 < yU 3 < —5 (2 flavors), —5 < /13 < —3 (6 flavors), 
—3 < fi-s < — 1 (6 flavors), — 1 < ^3 < 1 (2 flavors). We summarize it in Fig. [2j We note 
that nonzero d^ gives larger minimal-doubling range as shown above. 

What we want to study in this work is how these parameter ranges evolve in the finite 
gauge coupling direction. The most important question for practical use of this formulation 
in lattice QCD is how the two-flavor range, or the minimal-doubling range changes in the 
interacting theory. Our question is deeply related to a possible chiral phase structure with 
respect to the U(l) chiral symmetry. It is because spontaneous chiral symmetry breaking is 
expected to occur only in parameter ranges with physical fermions. We thus speculate that 
boundaries between x^SB and non-^SSB phases starts from fi% = —7 and ^3 — 1 in the 
weak-coupling limit. From next section, we will elucidate the fi3-g 2 chiral phase diagram by 
using strong-coupling lattice QCD and the Gross-Neveu model. 
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FIG. 2: The number of species in Karsten-Wilczek fermion as a function of the parameter On 
the boundaries between different sectors, the dispersion relation of fermions becomes unphysical. 

III. STRONG-COUPLING LATTICE QCD 

In this section we employ the strong coupling analysis to investigate the chiral phase 
structure in lattice QCD with Karsten-Wilczek (KW) fermion. The first step is to derive an 
effective potential of meson fields corresponding to the fermion action in Eq. ((8]). The strong- 
coupling study for KW fermion was first performed by the present author and collaborators 



m [38(. We here take the same approach. 

Lattice fermion action is generally written as following by using hopping operators P^- 
and an onsite operator M 

S = J2 MP^n + , ~ P-^n-v) + ^nMlpn. (12) 

By usingthese operators an effective action for mesons in the strong coupling limit can be 
written [8| as 



S eS (M) = N c ^ 



^2 Tr / ( A n>M ) + tr MM (n) - tr log M (n) 



(13) 



n, M ^ , V ) N 



c 



vi = ^p-c+a > V& = -^p;r n , (14) 

Tr/(A^) = -tr f (-M(n)(P+) T M(n + fi)(P-) T ) , (15) 
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where N c is the number of colors, Tr ( tr ) means a trace over color(spinor) index, and 
Ai(n) is a meson field. a,b are indices for colors while a,f3 for spinors. The explicit form 



of the function / is determined by performing a one-link integral o 
large N c limit, it is known that f(x) can be analytically evaluated 

1 + v / TT4x 



the gauge field. In the 



3 



as 



f{x) = VTTkc - 1 - In ^- = x + 0(x 2 ). (16) 

For most cases of studying phase structure, we can approximate it as f(x) ~ x, which 



corresponds to a large-dimension limit [39|. Since the phase transition is expected to be 
second-order for a massless case, this approximation at least works well near the phase 
boundary. We however note that it becomes less valid for large a. In the case of the 
Karsten-Wilczek fermion, we have M = 777.I4 + i(fis + 3r)7j and 

p+ = I |(7m + ir l4) /i = 1, 2, 3 p- = J - zr 7 4 ) H = 1, 2, 3 
I §74(1 + ^4) At = 4 I ±7 4 (l + d 4 ) At = 4 

We here assume a form of meson condensate with chiral and 4th vector condensates as 

M = crl 4 + OT474. (18) 

It is because the flavored chemical potential term is expected to produce 4th vector conden- 
sate, which is related to quark density. (We will discuss possibility of other condensates in 
the end of this section.) The explicit form of the effective action for a and 7r 4 is given by 

SW = -4Ar c Vol.Veff(a,7r 4 ), (19) 
V e ff(<7, 7r 4) = ^ log(cr 2 + vr 2 ) - ma + (/i 3 + 3r)7T 4 

- ^[3(1 + r 2 ) + (1 + eZ 4 ) V " ^[3(1 - r 2 ) - (1 + rf 4 ) 2 ]vr 2 . (20) 
We now find saddle points of S e s(A4) from 

— — = — — = 0. (21) 

0(7 07T4 

Then gap equations are given by 

3(l+r 2 ) + (l + d 4 ) 2 o 



a + m--^ = 0, 



2 a 2 + 7i 



3(l-r 2 )-(l + rf 4 ) 2 , , . . tt 4 n 

ff4 _(^ + 3r)--^-- ! = 0. (23) 
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It is notable that these gap equations have a particle-hole symmetry as (7r 4 ,/i 3 + 3) -H- 
(— 7r 4 , — yU 3 — 3), which is reflected by chiral phase structure as we will see later. We first 
consider m — 0, and solve the equations analytically. One of the main purposes here is to 
find a boundary between chiral symmetric and broken phases. For this purpose we take 
o" = after dividing the first equation by a since a is an order parameter of chiral symmetry 
breaking. Then we have 

3(l + r 2 ) + (l + d 4 ) 2 1 



2 ' 



(24) 



3(l-r 2 )-(l + d 4 ) 2 , _ 1 



vr 4 - (/i 3 + 3r) = — . (25) 

Z 7T 4 

These equations give chiral boundaries for /i 3 as 

6r 2 + 2(l + rf 4 ) 2 

u 3 = ± . = - 3r. (26) 

v/6r 2 + 2(1 + d A ) 2 + 6 

Therefore we have two ranges of /x 3 with different chiral properties, I and II : 

6r 2 + 2(l + d 4 ) 2 o 6r 2 + 2(l + d 4 ) 2 

I : /x 3 < v ' - 3r, /x 3 > v y - 3r, 27 

V6r 2 + 2(1 + rf 4 ) 2 + 6 v/6r 2 + 2(1 + d 4 ) 2 + 6 

II: ^ + 2(l + d4 y , 6, 2 + 2(l + 4) 2 3r 

v/6r 2 + 2(1 + rf 4 ) 2 + 6 v/6r 2 + 2(1 + d 4 ) 2 + 6 

The question is which corresponds to the chiral symmetric or broken phases. We for a 
while look into the r = 1 case to show details of chiral phase structure. Since the change of 
(i 4 just shift the chiral phase boundary ( 126]) . we for a while take — to catch the chiral 
phase structure although we note that we need to take care of it for a Lorentz symmetric 



continuum limit. For r = 1 and <i 4 = the boundaries are given by /2 3 = /x 3 + 3 = ± a/32/7 
(/i 3 ~ —5.14, —0.86). Here we defined shifted /i 3 as /i 3 = yU 3 + 3. By solving the gap equations 
in Eq. (l8l)( )23|) for r = 1 and <i 4 = we derive two solutions of chiral and 7r 4 condensates as 



Mi: a = 0, vr 4 = -fl 3 + ^/i 2 - 2, (29) 

and 

Atf : a = , 4 = _|. (30) 

By comparing effective potentials of the two solutions, we find which corresponds to vacua for 
the two parameter ranges, I and II. We show the following by substituting the two solutions: 



V eS (M£) - V eff (.M B ) > for -v/32/7 < fl 3 < a/32/7 (II) while V cG (M£) - V cS (M*) < 
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FIG. 3: Chiral condensate in strong-coupling lattice QCD with KW fermion for m = and r = 1. 
We define fis = fis + 3 and depict the chiral condensate as a function of £13 for < ^3 < 4. This 
result indicates the second-order chiral transition although the approximation f(x) ~ x becomes 
less reliable for larger a. 

for ^ 3 < -VW?, 

/U 3 > ^32/7 (I). To sum up, the chiral symmetric and broken phases 

are given by 



a = for /13 < -y/32/7 - 3, /x 3 > a/32/7 - 3 (I), (31) 



a = V 32 7{ ^ 3 + 3)2 for -v^277-3</i 3 < vW^-S (II). (32) 
4a/ 7 

The behavior of chiral condensate is depicted as a function of jx$ = 7x3 + 3 in Fig. [31 which 
shows that the transition between chiral symmetric and broken phases is second-order. We 
also depict parameter regions for the two phases in the strong-coupling limit for r = 1, 
c?4 = and m = in Fig. HJ As seen from Eq. (|26|) . a nonzero c/ 4 gives a larger physical 
range (II) with SSB of chiral symmetry. 

We show that the order of transition becomes a crossover for m 7^ 0. Figj5] shows the 
chiral condensate as a function of JI3 = [13 + 3 and m. It is obvious that the second-order 
phase transition changes into a crossover for m^O. We note that ir^ condensate is in general 
non-zero in the chiral-broken phase. 7T 4 becomes zero only for /x 3 + 3 = 0, which corresponds 
to a theory with unphysical dispersions. 
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<G>=0 



9 = °°- 



-3-V(32/7) 



<a>#) , <o)=0 
1 H -//3 



-3+V(32/7) 



FIG. 4: Chiral phase structure in the strong-coupling limit of lattice QCD with KW fermion for 
r = 1, di = and m = 0. 




FIG. 5: Chiral condensate in strong-coupling lattice QCD with KW fermion for r = 1 and g?4 = 
as a function of /J3 and m. 

We next look into mass of mesons within the same framework. We here consider m = 
but general values of r and d±. To calculate meson masses we expand the meson field as, 



M in) = Ml + * x (n)r T x , Xe{S, P, V a , A a , T Q/3 } 



(33) 



where S,P,V a , A a and T a/3 stand for scalar, pseudo-scalar, vector, axial- vector and tensor 
respectively. Here Ai is the vacuum expectation value of Ai(n). We note that 

1 4 



75 -p 7a -p ^757ct -p 

1 s - "2", ip - y, 1 v Q - -y , 1 a q - — 2" — , 1 r, 



a/3 



7a7/3 / ^ m 



(34) 
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Then the effective action at the second order of n x is given by 



N c I -0^K X (-p)D XY (p)7r Y (p), (35) 



where 



D XY (p) = ^(D XY (p) + D YX (-p)), (36) 

D XY (p) = itr(^ 1 r x >io 1 r y ) + ^tr(r x p / ;r y p / |)e^. (37) 



2 

In our case A4q = crl + £7^74 gives 



M 1 = o 1 2 ( ffl " Z7r 474). (38) 



<7 2 + 7r| 



We now write the whole inverse meson propagator matrix in the S-V-T-A-P sector as 
Dsvtap = 



t~\ s~i —irs\ —irs2 —irs3 



-c d Va 

irs2 n 
~2~ Ut ^ 



D 



-irs2 irsi 



V3 2 2 



D 



y i 2 2 

D Tl2 C 



u Tl3 — 
D T „ C =^sx 



'T2Z 

irs2 —irss /~i r-> 

— U U M 

—irsi irs3 r~i y\ 

—3— — -o U A2 

irsi —irs2 r~i 7~i 
— —3— t> i^A 3 

irs3 —irs2 irsi j~\ 
2 2 2 ^ 



PA 4 

(39) 
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1 



2(<X 2 + 7T|) 
1 

2((T 2 + 7T 2 ) 
1 

2(a 2 + tt|) 
1 

2(a 2 + 7r|) 
1 

2( C 7 2 + vr|) 



where components are given by 

D s 

D V4 
D Vs 

D V2 

D Vl 
D Tli 
Dr 2i 
Dt 34 
Dt 12 
Dt 13 
Dt 23 

D Al 

Da 2 

Da, 

D Ai 

D P 
C 



<J — tt: 



2(a 2 + 7T 2 ) 



2\2 



cr - 7r; 



2((T 2 + 7T 2 



2\2 



2\2 



2(a 2 + 7T 2 ) 

2 



2 2 



2(cr 2 + 7T 2 ) 2 
1 

2(a 2 + vr 2 ) 
1 

2(a 2 + 7T 2 ) 





-r 2 


)(ci + c 2 + c 3 ) + (1 + d 4 ) 2 c 4 ], 


(40) 




-r 2 


)(ci + c 2 + c 3 ) - (1 + d 4 ) 2 c 4 ], 


(41) 




r 2 ) 


;ci + c 2 ) + (r 2 - l)c3 + (1 + d 4 ) 2 c 4 ], 


(42) 




r 2 ) 


;ci + c 3 ) + (r 2 - l)c 2 + (1 + d 4 ) 2 c 4 ], 


(43) 




r 2 ) 


% + c 3 ) + (r 2 - l)ci + (1 + d 4 ) 2 c 4 ], 


(44) 


If/ 2 


-1) 


+ c 3 ) + (1 + r 2 ) Cl + (1 + d 4 ) 2 c 4 ], 


(45) 


--[(r 2 - 


-1) 


(ci + c 3 ) + (1 + r 2 )c 2 + (1 + d 4 ) 2 c 4 ], 


(46) 




-1) 


(ci + c 2 ) + (1 + r 2 )c 3 + (1 + rf 4 ) 2 c 4 ], 


(47) 


X r/ 2 
" + 4 [(r 


- 1 


)(ci + c 2 ) + (1 + r 2 )c 3 + (1 + rf 4 ) 2 c 4 ], 


(48) 


X r/ 2 
- + i [(r 2 


- 1 


)(ci + c 3 ) + (1 + r 2 )c 2 + (1 + 4) 2 c 4 ], 


(49) 




- 1 


)(c 2 + c 3 ) + (1 + r 2 )ci + (1 + d 4 ) 2 c 4 ], 


(50) 




-r 2 


)(c 2 + c 3 ) + (r 2 - l)ci + (1 + 4) 2 c 4 ], 


(51) 


lr/ 9- 
-+ 4 [(1+- 2 


)(ci + c 3 ) + (r 2 - l)c 2 + (1 + rf 4 ) 2 c 4 ], 


(52) 


lr/ 9- 
-+ 4 [(1+- 2 


)(ci + c 2 ) + (r 2 - l)c 3 + (1 + rf 4 ) 2 c 4 ], 


(53) 




r 2 ) 


;ci + c 2 + c 3 ) - (1 + c? 4 ) 2 c 4 ], 


(54) 




r 2 ) 


^ci + c 2 + c 3 ) + (1 + d 4 ) 2 c 4 ], 


(55) 
(56) 



(<x 2 + vr 2 ) 2 ' 

with Sk = sinpk and Ck = cospk- By diagonalizing this matrix, we can derive an explicit 
form of physical meson propagators. We note that C gets zero in the chiral symmetric phase 
while it has nonzero values in the chiral broken phase. 

We first check that pion mass becomes zero in the chiral-broken phase. For this purpose we 
substitute p = (0, 0, 0, imp) into the propagator. Then D P is decoupled and the calculation 
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gets simplified. We note that, in the chiral broken phase, the gap equation Eq. (jSJ) gives 
a 2 + 7I4 = 2/ [3(1 + r 2 ) + (1 + d 4 ) 2 ]. Pion mass in the chiral broken phase with r = 1 and 
m = is derived from a pole of Dp with p = (0, 0, 0, imp) as 

£> P (0, 0, 0, imp) = { — LZ^L^l - + d 4 ) 2 cos(mp) + 3(1 + r 2 )) = 

— > cosh (mp) = 1. (57) 

This mode corresponds to a massless NG boson associated with SSB of chiral U(l) symmetry, 
which is consistent with realistic QCD. 

On the second-order phase boundary, a divergent correlation length in the critical behav- 
ior should produce another massless mode. We here show that the scalar meson becomes 
massless on the critical line. An inverse meson propagator matrix in general has off-diagonal 
components, but substitution of p = (0,0,0,ims) again decouples the scalar sector in the 
chiral-symmetric phase (C = 0). An explicit form of the scalar inverse propagator is given by 
Eq. (140]) . Then we derive the scalar mass from a pole of this propagator with p = (0, 0, 0, im s ) 
as 



£> s (0,0,0,im s ) = ->• coshm s = l- 1 + 7 f— — L + \ *' , (58) 



3(1 + r 2 ) 1 2(a 2 - ir 2 

(l + d 4 ) 2 + (l + d 4 ) 2 (a 2 + 7r 2 ) 2 



As seen from this, scalar meson has nonzero mass in general. On the phase boundary, how- 
ever, we have a = 0, 7r 4 = — u 3 + + 2[3(1 — r 2 ) — (1 + dA 2 ] and u 3 = ± , 6r +2( - 1+d ^ 

,4 P'A v "3 I \ J ^ v /6r2+2(l+d 4 ) 2 +6 

with fi 3 = fi 3 + 3r. It leads to zero a mass as 

cosh(ms) = 1. (59) 

We can also show this from the chiral-broken phase. We have two massless modes only on 
the 2nd-order phase boundary, a-meson and 7r-meson, which is inconsistent with the meson 
mass spectrum in QCD. We thus need to avoid this point in the simulation of QCD. In the 
first place this boundary corresponds to the boundary between the two-flavor and no-flavor 
ranges, which has the unphysical dispersion relation as shown in Sec. [Til From the theoretical 
viewpoint, it is however an attractive topic. Further study including numerical simulations 
can elucidate detailed properties of this point. 

Now let us discuss possibility of restoration of Lorentz symmetry. In the strong-coupling 
limit, what we can do is tuning of [13 and g? 4 since we can ignore tuning for the plaquette 
action in this limit. However, we cannot restore the Lorentz symmetry correctly in this 
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limit because of the large lattice artifacts. The A 4 meson propagator explicitly manifests 
this point as following: The v4 4 -sector propagator, which is diagonal even with finite spacial 
momentum, is given by Eq. ( 15^1) . In the chiral-broken phase, we have l/(cx 2 + 7r|) = [3(1 + 
r 2 ) + (1 + gLi) 2 ]/2. By substituting p = (0, 0, 0, irriA^), the pole of the propagator determines 
the axial vector meson mass as 

coshm ^ = 1 + (IT4F (60) 

For small but finite spatial momentum as p = (pi,p 2 ,p3,iE), the pole of the propagator up 
to 0(p 2 ) gives 



r 



2 



B2 = < + (TW p2 ' (61) 

It apparently indicates that we need to tune as (1 + d^) 2 = 1 — r 2 . However, for r = 1, 
the dependence on p itself disappears and the Lorentz symmetry cannot be restored. It 
is a typical strong-coupling artifact. Moreover, the tuned value of d^ does not depend on 
/i3, but it is also a lattice artifact. In the first place, the equation flB]) indicates that tt^ 
always has a nonzero value in the physical parameter range (a ^ 0) even if we tune /13 
and di independently except for fi^ + 3r = 0. All these results show that it is difficult to 
restore the symmetry in the physical parameter range within the framework of the strong- 
coupling QCD. We expect that it is just a strong-coupling artifact, and we can make ir^ 
zero and restore Lorentz symmetry in the weak coupling by tuning the three parameters 
appropriately. 

In this section we have assumed the form of condensation a + 274^4. As shown in the 
Appendix. A, we can also consider possibility of other condensations as o + 274^4 + 275^5 or 
a + Z747T4 + ^7475^45. The results show that the solution with nonzero or 7145 condensates 
cannot be a vacuum, and our solutions in this section are likely to be true vacua. We 
thus consider that the parity breaking phase will not appear in the KW fermion unless 
we introduce the flavored-mass terms shown in 32|. At the weak coupling there may be 



a more subtle competition between the discretization error and the counterterm. We thus 
need further study to conclude whether the parity breaking exist or not when we take into 
account all the three counterterms in the weak coupling. 

In the end of this section, we discuss the other type of minimal-doubling fermions, called 



the Creutz-Borici type 15|, |l6[. We can analyze it in a parallel way. We note that this 



type specifies the diagonal direction characterized by 2r = 71 + 72 + 73 + 74, instead of 
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the time direction. Appendix. B is devoted to detailed analysis for this case. The result 
is qualitatively the same. We find two chiral boundaries and chiral condensate is nonzero 
between the boundaries. 



IV. GROSS-NEVEU MODEL 

We investigate the whole phase diagram for Karsten-Wilczek (KW) fermion by using the 



two-dimensional lattice Gross- Neveu model 40M49| . which has common features with 4d 



lattice QCD. In two dimensions, massless KW fermion action is given by 

2 

aD K w(p) = i 7„ sin ap u + 272 [r(l — cosapi) + ji + d sin ap 2 }. (62) 

In this section we concentrate on the case with r = 1. In this section we denote the relevant 
parameter /i 3 as just fi = fi 3 and denote the maginal parameter d± as d = d±. To look into 
the number of flavors in a free theory, we for a while consider d — 0. For —3 < \i < — 1 and 
— 1 < \i < 1, there are only two zeros, and it becomes minimal-doubling. For \i < —3 and 
/j, > 1, there is no zero, and it becomes a fermion- less theory. For \x = —3, 1, there is one zero, 
but the dispersion relation becomes unphysical ~ p\ + p\. For fi = —1, there two zeros but 
whose dispersion relation is again unphysical. The main difference from four-dimensional 
cases is that there is no 6-flavor range. 

The lattice Gross-Neveu model with KW fermion is given by 



n,v n 
n 

+ ^i~f2*Pn + VV72(^ n+2 - - + m £ Mn, (63) 

n n n 

where v stands for v = 1,2, n = (^1,^2) are the two dimensional coordinates and ip n 
stands for a iV-component Dirac fermion field {ip n )j{j = 1, 2, ...,N). The bilinear ijjijj means 
YljLi i^j^ii an< i {da-i 9%) corresponds to the 't Hooft couplings for the two types of four-fermi 
interactions. Here we define the two dimensional gamma matrices as 71 = a±, 72 = 02 and 
73 = 03. We make all the quantities dimensionless in this equation. We consider scalar and 
time-direction vector four-fermi interactions, which are natural choices since KW fermion 
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y^[^(u(n) - m) 2 + \(n 2 (n) - /i - l) 2 ] + ^^„[cr(ra) + i^ 2 iT 2 (n)]tp n 



specifies the time direction. We note that, if we drop a mass term (m = 0), the action has 
Ti / 2 discrete chiral symmetry, which can be spontaneously broken due to chiral condensate. 
By introducing auxiliary bosonic fields a(n), ^(n) we remove the four-point interactions as 

n,u n n 

Z n 9a 92 

(64) 

By solving the equations of motion, we show the following relation between these auxiliary 
fields and the bilinears of the fermion fields 

a(n) =m- ^ip, (65) 
7r 2 ( n ) = 1 + n - ttLfay^ (66) 

These relations indicate that a and ir 2 stand for the scalar and vector mesons. After inte- 
grating the fermion fields, the partition function and the effective action with these auxiliary 
fields are given by 

Z = f Y[da(n)d7i 2 (n)e-* rs <*(°™\ (67) 
S eS (a,n 2 ) = ^[4(^(n)-m) 2 + ^(7r 2 (n)-^-l) 2 ]-Trlo gj D n , m) (68) 



it 



9 a 92 



with 



D n , m = [cr(n) + 172^2 W]fl n .m + — (5„,+ M , m - 5. 



2 \-n+[i,m u n—fi,m) 
il2 (K , r \ , ,72 

Y^ d n+i,m + d n-i,m) + d ^ 



-TriK+hm + S n-lm) + rf ir( 5 n+2,m ~ S n-2,m)- ( 69 ) 



Here Tr stands for the trace both for the position and spinor spaces. As is well-known, the 
partition function in the Gross-Neveu model is given by the saddle point of this effective 
action in the large N limit. We denote as <j(n), n 2 (n) solutions satisfying the saddle-point 
conditions 

5S cS [a(n),7r 2 (n)] SS eS [a{n), 7T 2 {n)] = 

Sa(n) 5n 2 (n) 1 ' 

Then the partition function is given by 



-NS c ff [a,if2] 



(71) 
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By assuming the translational invariance we define the position-independent solutions as 
a = cr(0) and 7r = 7f 2 (0) Then we can factorize a volume factor V = ^ n 1 in the effective 
action as 

S cS = VS cS (a,n 2 ), (72) 

Ses(<T,ir 2 ) = =-5(0 - mf + -^(7r 2 - fi - l) 2 - -^Tr log D. (73) 
^9(7 ^92 v 

We can write Tr log D in a simple form by the Fourier transformation to momentum space 

/d 2 k 
^log[det(cr + ry 2 7r 2 + ry 2 ((l + d) sin k 2 - cos k x ) + sin fci)] 

/d 2 k 
^plog[a 2 + (vr 2 + (1 + d) sink 2 - cosh) 2 + (sin fci) 2 ], (74) 

with det being the determinant in the spinor space. Now saddle-point equations are written 



as 

f_ 

5a gl J (27r) 2 a 2 + (n 2 + (1 + d) sin k 2 — cos ki) 2 + (sin k\) 2 

5S e ff n 2 — ji — 1 f d 2 k 7r 2 + (1 + d) sin k 2 — cos k\ 



^ {a ~ m) -*I^ ^TTZ ,,: J ^ , v , /: , , ,, 0. (75) 

0. (76) 



5tt 2 g 2 J (27r) 2 ct 2 + (tt 2 + (1 + d) sin A; 2 — cos &i) 2 + (sin k\) 2 

Here the values of a and 7r 2 in the vacuum are determined as a(fi,d, g 2 , g 2 ,m), 
7r 2 (/i, d, g 2 ., g 2 , m) from the saddle-point equations once /i, d, g 2 , g\ and m are fixed. 

Let us look into the phase diagram with respect to chiral symmetry. We here consider 
a massless case as m = to have the exact discrete chiral symmetry in the action. To 
capture rough structure of the phase diagram, we first take the simplest case with d = 
and g 2 = g 2 = g 2 . Nonzero d just gives slight change of the phase diagram. Since a single 



coupling constant works when we study the Aoki phase in the Wilson fermion [27J , we expect 
that the above condition for the couplings works at least for deriving rough phase structure. 
The order parameter is a, which can be zero or non-zero depending on values of [i and g 2 . 
The phase boundary is determined by imposing a = on Eq. fl75|)fl76|) after the overall a 
being removed in Eq. (1751) . Then the conditions for the phase boundary are given by gap 
equations as 

7r 2 — u c — 1 f d 2 k ti 2 + sin k 2 — cos ki , , 

(77) 



g 2 J (27r) 2 (7r 2 + sin k 2 — cos ki) 2 + (sin ki) 2 ' 

g 2 J (27r) 2 (7r 2 + sinfc 2 -cosA; 1 ) 2 + (sinA; 1 ) 2 ' 1 ' 
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FIG. 6: Chiral phase diagram for KW Gross-Neveu model in the fi-g 2 space. A stands for the 
chiral symmetric phase. B stands for the chiral-broken phase. The chiral boundaries are connected 
to the edges of minimal-doubling ranges. 

with fi c being a critical value of \x. Here we derive the chiral phase boundary fi c {g 2 ) as a 
function of the coupling g 2 by getting rid of vector condensate tt 2 from these equations. The 
phase diagram is depicted in Fig. [6j A stands for the chiral symmetric phase a = and 
B for chiral broken phase a ^ 0. For — 2 < 7r 2 < 2 the integral in Eq. (l78j) diverges, which 
leads to the critical lines with zero gauge coupling for —3 < fi < 1 as shown in Fig. |6j 
This is reasonable since the weak-coupling limit should have zero chiral condensate. As we 
expected, the chiral critical line is connected to boundaries between two-flavor and no-flavor 
phases in the weak-coupling limit (fi = —3, 1). It is consistent with our intuition that the 
fermion-less theory cannot cause spontaneous breaking of chiral symmetry. We can check 
the mass of the scalar meson mass becomes zero on the critical line fi c (g 2 ). The mass of o 
is calculated analytically on the boundary, and is shown to be zero as 



5a(n)5a(m) ,flc 5 2 a 2 lAt 
V 



i r d 2 k i 

— — — 2 



g 2 J (2tt) 2 a 2 + (7r 2 + sin k 2 — cos ki) 2 + (sin k\) 2 
d 2 k 1 



-4a 2 



(27r) 2 a 2 + (ir 2 + sin k 2 — cos ki) 2 + (sin fci) 
0. (79) 
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A massless scalar meson indicates that the phase boundary we derived is a second-order 
critical line. This result is consistent with the strong-coupling lattice QCD in the previous 
section. 

Next, we discuss more general cases with d 7^ and g 2 7^ g\. Nonzero values of d change 
sinfc 2 — > (1 + d)sink 2 in (1771) (1781) and give just qualitative changes of the phase diagram: 
Since the minimal-doubling ranges are given by — 1 < fi < 1 + d and —3 — d < \i < — 1 for 
nonzero d, it gives a larger physical range in the phase diagram. As an example, we depict 
the n~g 2 phase diagram for d = 0.5 in Fig. [7J In the case with two independent coupling 
constants g 2 7^ g 2 , the equations for the phase boundary are given by 

, _ A g%\ _ 2 f tt 2 + sin k 2 - cos k x 

^ + 2 V 9 2 J 92 J (2^) 2 F2 + sin k 2 - cos k,f + (sin k,) 2 ' m 

1 _ o [ d 2 k 1 

gl J (27r) 2 (7r 2 + smk 2 -cosk 1 ) 2 + (smk 1 ) 2 

In this case the phase diagram is deformed to some extent: If we fix one of the coupling 
constants and depict the phase diagram, the physical phase gets larger toward the strong- 
coupling limit while it gets narrower at the weak-coupling. As an example, we depict the 
fi-g 2 phase diagram for g\ = 1.0 in Fig. |HJ Since we have already shown that the physical 
phase gets rather narrower in the strong-coupling limit in 4d lattice QCD in Sec. IHII it is 
natural to consider that the single- coupling GN model in Fig. [6] is sufficient to mimic 4d 
QCD at least for investigating the chiral phase diagram. We note that we can introduce two 
more types of four-point interactions in the GN action (1631) as (^n^TiV'n) 2 and (^ n ^73'0n) 2 \ 
which give tt\ and ir meson fields. Both vacuum expectation values of tti and 7r, if they are 
nonzero, cause spontaneous parity symmetry breaking. However, by comparing the effective 
potentials or solving four gap equations (^jf- = 0, = 0, = 0, ^jf- = 0), we can 
clearly show that they have zero vacuum expectation values as in the case of the strong- 
coupling lattice QCD in Appendix. A. Instead of solving all the gap equations numerically, 
we here show the consistency check on 7Ti = 0. From two saddle-point equations = 
and -j^- = we derive one equation, 



d 2 k sin fci 
q _ 2 / (82) 

(27r) 2 a 2 + 7r 2 + (ir 2 + (1 + d) sin k 2 - cos ki) 2 + {-k x + sin ki) 2 ' 



If we take 7Ti = 0, the right-hand side becomes an odd function of k\, and the above equation 
holds identically. We can also show from the gap equations that it = is preferred. These 
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FIG. 7: Chiral phase diagram for KW Gross-Neveu model in the [i-g 2 space with d = 0.5 (g 2 = 
g 2 = g 2 ). The physical (minimal-doubling) phase with chiral condensate is enlarged by nonzero 
values of din the weak-coupling limit, the minimal-doubling phases are given by —3.5 < fi < — 1 
and — 1 < fx < 1.5. 

results mean that the parity-broken vacuum is not preferred, and as long as we consider the 
vacuum of the theory, the two four-fermi interactions in fl63|) are sufficient. We next show 
that, however, we have to take into account the excitation of 7Ti and it to investigate the 
Lorentz symmetry restoration in the GN model. 

By now we have considered only the vacuum of the theory. Although it works to elucidate 
the chiral phase structure, we need look into excitations from the vacuum in order to in- 
vestigate renormalization for the rotation symmetry restoration. By expanding the mesonic 
action up to the 2nd order, we can derive the meson propagator matrix for a, 7Ti, vr2, 7t as in 
the strong-coupling lattice QCD in Sec. Ill IL By diagonalizing this 4x4 matrix, we obtain the 
proper propagators. Then, we can discuss how the rotation symmetry can be restored by 
tuning the parameters including /x, d, g 2 and g\. Although we can perform this procedure in 
principle, all the matrix components include complicated momentum integrals in this case. 
Moreover, we need to substitute the VEV of meson fields derived from the gap equations, 
which also require numerical integrals. We consider that these numerical calculations are 
beyond the scope of this study, but the future work should be devoted to it. Instead, we here 
show a process of deriving the dispersion relations in details in the following. The meson 
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FIG. 8: Chiral phase diagram for KW Gross-Neveu model in the [i-g^ space with g\ = 1.0 (d = 
0). The physical (minimal-doubling) phase with chiral condensate gets larger toward the strong- 
coupling limit. 

excitations from the vacuum are given by 

o = a + a(n), m = 7fi + 7Ti(n), 7r 2 = 7f 2 + T^fa), n = Tx + n{n) 1 (83) 

where VEVs (a, Tti, 7T 2 , vf) are determined by the gap equations. We have already shown 
7fi = 7f = while a and 7t2 depends on the parameters. The Dirac operator with these 
excitations are written as 



n,in 

+ M(n)S n 
i T 

(<Sq )n,m = [&+ il2^2\5 n .m + -£(Sn+n,m ~ &n-ii,m) / & 

M(n) = a(n) + z'7i7Ti(n) + i^2(n) + 27 3 7r(n), 



(84) 

(85) 
(86) 



where we make the lattice spacing a manifest. Note a, 7Ti, 7T2, 7t ~ 0(1/ a). For completeness 
we consider the following action including four independent coupling constants, 

n \ N r / \2 2 ' - ^ ° 



— 7r l 2 + 7772( 7r 2-^-l) 2 + ^7T 2 ]-Tr log[5 - 1 (l+5 7W)]n,m. (87) 



Tr is a trace for coordinate and spinor spaces. By expanding this expression up to second 
order of mesonic fluctuations, we derive the effective action up to the quadratic order as 



C _ q(0) , o(2) 



(88) 
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with 



S. 



(2) 
off 



E 



a(n) 2 7Ti(n) 2 vr 2 (n) 2 ir 2 (n) 1 v-^ „ , s . ^ , s . 
-^ J - + + + ~yy- + rtr Sb(m, n)A4(n)S„(n, m)A<t( 



2^ 



2<7? 



2^ 



(89) 



with tr is a trace for the spinor space. By fourier transformation we can derive the form in 
the momentum space as 



S, 



(2) 
eft' 



n/a d 2 p 



K/a (27T) 2 L 2gl 2g 2 <lg\ 



~tr ^S (p+k)M(p)S (k)M(-p) 
2 A/a (2tt) z 



with 



<So(k) 



det(<S - 1 (A;)) = cf 2 + 



O" - ?72 7T2 



sin a/c 2 cos a/ci 
a) 1 



7T 2 



d) 



sin ak 2 cos a/ci 



a a 

Finally, the mesonic propagator matrix is given by 



+ 



sin ak\ 



271- 
2 



sin a/ci 



V X y{p) 



2c(2) 



5 2 S 



eft 



(91) 
(92) 

(93) 



&Mx(p)£My(-p) 

where X and Y stand for one of channels a, ttx, 7r 2 and it. For example, T> acj (p) is given by 

m , , i r /a ^ i 



& A/a (27iydet[S^(p + k)S^(k)} 



a 2 



- (vr 2 + (1 + d) 



sin a/c 2 cos ak\ 



vr 2 



d) 



sin a(p 2 + fc 2 ) cos a(pi + k\) 



sin afci sin a(pi + &i 



(94) 



In the same way we can derive all other 16 components of the matrix, all of which take 
nonzero values in general. What we are interested in is terms of the components up to 
0(p 2 ) as far as we consider the rotation symmetry up to 0(a) discretization errors. The 
coefficients can be extracted through numerical integrals of the equations as (|94|) . We note 
that we also need to substitute into the integral the values of VEV of o and 7r 2 derived from 
the gap equations. Then we diagonalize the 4x4 meson matrix T>xy to derive the proper 
meson propagator. By introducing p = (p, iE) into the propagators, we find the dispersion 
relations. The question is which parameters among /i, d, g 2 , gf, g 2 , g 2 need to be tuned 
to recover Lorentz symmetry. It can give an important suggestion to lattice QCD with the 
minimal-doubling fermions. We devote a future work to this analysis. 
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FIG. 9: Conjecture on fi^-g 2 chiral phase structure for Karsten-Wilczek fermion with r=l. The 
width of the minimal-doubling range determines how hard it is to tune ^3. 

V. CONJECTURE ON PHASE STRUCTURE IN QCD 

From the study of strong-coupling lattice QCD and the Gross-Neveu model we speculate 
on the whole chiral phase structure in lattice QCD with Karsten-Wilczek fermion. Fig. |9]is 
a conjectured chiral phase structure with the number of flavors in the [i^-g 2 space for r = 1. 
There are roughly two phases with and without chiral condensate, or equivalently with 
and without SSB of chiral symmetry. As was shown in the previous section, the boundary 
between chiral symmetric and broken phases starts from the edge of the two-flavor region of 
the free theory. We expect that the chiral boundaries are connected to the two-flavor and 
no-flavor phases also in 4d QCD as shown in Fig. [9j 

The question is a boundary between two-flavor and six-flavor ranges. In the weak-coupling 
limit (g 2 = 0) we analytically know the number of physical flavors: There are four sectors 
with two, six, six and two flavors. There are only no flavors of fermions outside these ranges. 
Toward the strong coupling, these ranges will change with g 2 as shown in Fig. [9j We have 
seen that we cannot distinguish two-flavor and six flavor ranges in the strong-coupling limit, 
which means that the number of species becomes an ambiguous notion in this limit. We 
thus expect that the boundary disappears at a certain gauge coupling, and the two-flavor 
and six-flavor regions become undistinguishable as shown in Fig. [9j 

From the viewpoints of practical application to two-flavor QCD, the relevant parameter 
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/i 3 has to be tuned to cancel the 0(1 /a) imaginary chemical potential renormalization for 
the two flavors. One necessary condition (but not a sufficient condition) for this purpose is 
to set in the minimal- doubling range to realize the two-flavor QCD. As we conjectured in 
Fig. [9j the minimal-doubling range in the middle gauge coupling should have some width. 
One possible indicator of the minimal-doubling phase is the pion spectrum: If /13 is in the 
no-flavor range, there is no SSB of chiral symmetry and no massless pion. If /x 3 gets into 
the six-flavor region, the number of pseudo Nambu-Goldstone bosons increases. However, 
setting /X3 in the minimal-doubling range is not sufficient for physical QCD to be described: 
The Lorentz symmetric dispersion relation is broken down for general values of ^3 in the 
minimal-coupling range as shown in Eq. (11 II) and below for a free theory. 



D(p) ~ ijiPi + (1 + di) 2 -fi 2 3 + 0(ap 2 ). (95) 

Since this free-theory argument indicates that the rotation symmetry can be restored by 
tuning g?4 as (1 + d^) 2 — fi\ — 1, we may be able to restore Lorentz symmetry just by tuning 
d A with fj, 3 being set in the minimal-doubling range. Note that the minimal- doubling range 
gets larger with nonzero d A as shown in Sec. Ill and IV, thus it seems that proper tuning of 
g?4 can be done for any value of /x 3 without breaking down minimal-doubling. (We need one 
more parameter tuning for the plaquette action in any case.) We also emphasize that the 
same relative tuning of ^3 and d A makes tree level couplings of the gauge field to the fermions 
have a correct Lorentz-symmetric form. To show this, we look into the quark-quark-gluon 
vertex at the tree level. For the case of [i^ = and d A = it is given by 

V(p, k) = -ig I 7 M cos ^ — z- + 74(1 - 5 M ) sin ^ — J , (96) 



as shown in 



For nonzero fi 3 and d^, it is modified as 




a( Pj + kj, 



V(p, k) = -ig ) 7j cos h 74 



a(p A + k A ) a (Pi + k i) 
sin 



1 + d 4 ) cos h 2 / 



2 

(97) 

where we have no direct emergence of /X3 since it is a parameter for onsite (non-hopping) 
terms. However, as we discussed, the zeros of the Dirac operator for nonzero /i 3 and d A is 
given by a function of /i 3 and d A as p = k = (0, 0, 0, - arcsin(— j^-))- Now we expand both 
p and k about the zeros as p — > p + p and k — >■ k + k. In particular, the coefficient of 74 in 
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97]) is expanded as 



a(p 4 + k 4 ) a(p 4 + k 4 ) . a(p 4 + k 4 ) . a(p 4 + k 4 )i v-^ . a(pi + kA 

cos cos sm sm + > sm 

2 2 2 2 



i=l 



'l + d 4 



3 

At3 a(p 4 + /c 4 ) /i 3 a(p 4 + /c 4 )l . a(pi + ^) 
cos h sm + > sm 



1 + d 4 ) 2 2 1 + d 4 



■i=i 



l + d 4 )Jl-—^— + 0(ap,ak), (98) 



[l + d 4 f 

Then, the vertex surviving in the naive continuum limit is given by 

V(p, k) = -ig ^71 + 72 + 73 + 741/ (1 + d 4 ) 2 - fi^j + 0(ap, ak). 



(99) 



Here we omit the ± sign in front of 7 4 for the doubler pairs for simplicity. It is now obvious 
that the tuning condition for the speed of light (1 + d 4 ) 2 — /i| = 1 also fixes the couplings of 
the gauge fields to the fermion fields in the tree level up to the discretization errors. At least 
in the naive continuum limit, we can have a correct set of the Feynman rules for fermion 
fields with the condition. However it is too early to conclude that this condition is sufficient 
for Lorentz symmetry restoration since all the other Ward identities may not be corrected 
by it in the interacting theory including the loop effects: To discuss details, we consider the 



quark self-energy in lattice QCD with minimal-doubling fermions following 26[ as 



74 

E(p,m) = ijuP^iip) + mH 2 {p) + ^i(fl'o) ■ hm + d 2 (g ) • i — , (100) 

a 

where Si, £2, d\ and d 2 can be calculated in the perturbative analysis. It is clear that /X3 and 
d 4 corresponds to counter parameters for d 2 and d\ respectively, and the Dirac operator fl95|) 
is renormalized as d 4 — > d 4 + d\ and fi^ — > fi^ + d 2 in the interacting theory. The last term 
with d 2 , or the 0(1 /a) renormalization, causes a shift of the poles of the Dirac propagator 
away from their original positions as well as the change of the speed of light as shown in 
( 195]) . In the present work, we have also shown that this contribution changes the size of 
the minimal-doubling range at the finite gauge coupling and provides the non-trivial phase 
structure as shown in Fig. HJ The question is whether or not we need to move the poles of 
the propagator back to the tree-level positions by tuning /i 3 for a correct continuum limit. 
As far as the dispersion relation can be restored by d 4 , it seems that the position of poles 
is not relevant to physics. However, in practical use of the minimal-doubling fermion, d 4 



should be also tuned to make the conserved charge unity as shown in 



261 ] . It is not obvious 
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whether this condition can also fix the speed of light non-perturbatively. More generally 
speaking, it is very nontrivial whether all the Ward identities are fixed only by one tuning 
condition in the interacting theory. If we cannot restore Lorentz symmetry only with 
tuning unlike the free theory, it means that we still need to fine-tune the three parameters 
independently for a correct continuum limit of lattice QCD simulations 24k26| . Further 
study is needed to figure out this point. 

In the end of this section, we comment on another possibility for studying minimal- 
doubling fermions. One interesting possibility is the chiral perturbation theory for minimal- 
doubling fermions. Although we expect that it is quite tedious to construct the minimal- 
doubling ChPT with the lower discrete symmetry than Wilson and stag gere d fermions, the 
process could have some similarities with that of the in- medium ChPT 50[. If we succeed 
to construct the minimal-doubling ChPT, it is intriguing to consider the Lorentz symmetry 
restoration within the theory and discuss the parameter tuning for the symmetry restoration. 
We can also investigate the vacuum and the phase structure in the theory. We devote future 
works to the study on the minimal-doubling ChPT. 



VI. SUMMARY AND DISCUSSION 



In this paper we investigate the chiral phase structure in the parameter space for lattice 
QCD with minimal-doubling fermions, which can be seen as lattice fermions with a species- 
dependent naive chemical potential term. We study the phase structure with Karsten- 
Wilczek fermion by using strong-coupling lattice QCD and the Gross-Neveu model, and find 
out the nontrivial chiral phase structure in the /X3 — g 2 plane. 

In Sec. mi we have proposed flavored-chemical-potential lattice fermions, where some of 
doublers are eliminated by a species-dependent chemical potential term without losing all 
chiral symmetries. Minimal-doubling fermions are shown to be a special case of this type. 
In Sec. IIHI we investigate the chiral phase structure of lattice QCD with Karsten-Wilczek 
fermions in the strong-coupling regime. We derive an effective action for the scalar and 4th 
vector fields, and find that chiral symmetry is spontaneously broken in a certain range of the 
relevant parameter /i3 while the chiral condensate becomes zero outside the range. We show 
that there is a 2nd-order phase transition between chiral symmetric and broken phases as a 
function of fi%. We also show that pion becomes massless as a Nambu-Goldstone boson in 
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the chiral-broken phase while the scalar meson becomes massless only on the second-order 
phase boundary due to the critical behavior. In Sec. II VI we obtain information on the whole 
chiral phase structure in the — space by using the Gross-Neveu model with KW fermion 
in large N limit. From the gap equations we derive a chiral phase diagram. In Sec IVl we 
discuss the whole phase structure in 4-dimensional lattice QCD with KW fermions. We 
conjecture the chiral phase structure and numbers of massless flavors in the phase diagram. 
We also discuss the fine-tuning process of /i^ in lattice QCD from the viewpoint of the 
minimal-doubling range of the conjectured phase diagram. 

In this paper we have investigated whether or not chiral symmetry is spontaneously 
broken depending on the parameters. Unlike Wilson parity-flavor breaking phase with the 
width m ~ 0(a 3 ) j^J, it seems that we do not have fine phase structures in the weak- 
coupling regime as shown in Fig. [9j It is because the chiral symmetry breaking is physical, 
and there is no competition between the relevant parameter and the lattice artifact, which 



did cause the Aoki phase or the Creutz-Sharpe-Singleton phase in Wilson fermion [29|, [30] . 
Although our analysis on chiral symmetry and that on parity symmetry in Wilson have 
some common features, they are on the different levels. 

In this work, we do not take much care of 7r 4 condensate nor its physical implication. As 
we discussed in Sec. Hill this condensate is likely to be related to the (pseudo) quark density. 
On the other hand, Eq. (jSJ) in the strong-coupling limit shows that 7r 4 is nonzero for any 
values of [13 and except for the unphysical point /13 + 3r = 0. We consider that it is 
just a strong-coupling artifact, and the condensate can be eliminated by fine-tuning all 
the three parameters appropriately in the weak coupling. There is possibility that another 
way of introducing the dimension-4 counterterm as an 0(1) imaginary chemical potential 
term (^^+4 — e'^ipn-i) m &y work since the flavored-chemical-potential term could 
generate 0(1) effective imaginary chemical potential too. 

In the end of this paper, we refer to a future work. We can apply the flavored-chemical- 
potential (FCP) fermions to finite-temperature and finite-density lattice QCD, and obtain 
finite-(T,/i) QCD phase diagram. In our next work 5l| . we will propose a new method to 



study lattice QCD in medium by using FCP fermion formulations. 
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Appendix A: Possibility of pion condensation 

We here assume a form of meson condensate with chiral, 4th vector and pion condensates 

as 

Mo = &I4 + 27T474 + 27T 5 7 5 . (Al) 

The effective potential for a, 7r 4 and ir 5 is given by 
V eS (cr, vr 4 , vr 5 ) = - log(er 2 + 7r 4 + nf) -ma+ (/i 3 + 3r)vr 4 

- |[3(1 + r 2 ) + (1 + eZ 4 )V - ±[3(1 - r 2 ) - (1 + d^nj - ±[3(1 + r 2 ) + (1 + d 4 ) 2 ]vr 2 . 



(A2) 



We now find saddle points of V e ff given by 

3(l + r 2 ) + (l + d 4 ) 2 



3(1 - r 2 ) - (1 + <P 2 



2 



a + m 



v"4 - (a*3 + 3r) - 



2 

3(l + r 2 ) + (l + cV 2 





a 




= o, 


(A3) 




h7T 2 


+ 7T 2 




7T 4 




= o, 


(A4) 


a 2 - 


h7T 2 


+ 7T 2 








= 0. 


(A5) 


a 2 - 


h 7T| 


+ 7T 2 



2 

We now consider a case for r = 1, <i 4 = and m = 0. We have two types of solutions as 7r 5 = 
and 7r 5 7^ 0. For 7r 5 = 0, we have the two solutions A4q and A4q in ()29l)(l30!) as we discussed 
in Sec. HH For tt 5 ^ 0, {53} § ive ^ 2 + tt 4 + vr| = 2/7 and vr 4 = -/2 3 /4 = -(/i 3 + 3)/4. 
By substituting them into ( 1A2j) the effective potential as a function of /13 is given by 

Vc fr (7r 5 ^0) = ^log^-i-^ 2 . (A6) 

On the other hand, for example, the effective potential for one of 7T5 = solutions M.q is 
given by 

V cff (7r 5 = 0,M%) = ±log^ - ± - (A7) 
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It is obvious V c ff (7r 5 ^ 0) ^ V c ff (7r 5 = 0, M.q) for any value of p, 3 . It indicates that the 7r 5 ^ 
solution is unlikely to be a vacuum. We thus conclude that, at least in this framework, there 
is no pion condensate or no spontaneous parity symmetry breaking. We also perform the 
same analysis for the condensate form (7I4 + 2^474 + 2^45745, and will find that ^45 7^ is 
not a vacuum of the theory. 



Appendix B: Creutz-Borici case 



A free action of Creutz-Borici fermion is given by 



SBC = ^2i\Y1 ^n(lpn+n ~ ^n-fi) + ^ ^ ~ 7/0 ( 2 ^« _ ^n+£ ~ ^n-/0 

n 



+ ic 3 ^j n Tip n + mi(j n ip n ) , (Bl) 

T = \j2^ 7; = r 7 ,r = r - 7 „ (B2) 

where C3 corresponds to ^3 in Karsten-Wilczek fermion and T satisfies T 2 = 1 and {T, 7^} = 
1. In this case we have onsite operator and projection operators as M = ml 4 + i{c% — 2r)T T 
and 

P+ = i{ 7 „(l + ir) + irT}, P" = I{ 7 „(1 - ir) - irT}. (B3) 

The strong-coupling analysis is done in a parallel manner to KW fermion. By taking M. = 
a + Z7rrr, the corresponding gap equations become 

2(1 + r 2 )a + m- - - 2 = , (B4) 
cr 2 + 7Tf 

(1 + r 2 )vr r - (c 3 + 2r) - ^3 = . ( B5 ) 
For m = 0, equations for second-order phase boundaries are given by 

2(l + r 2 )-^ = 0, (B6) 

7T r 

(l + r 2 )vr r -( C3 + 2r)-— = 0, (B7) 

7r r 



Then chiral boundaries are given by 
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For r = 1 chiral condensate is nonzero and chiral symmetry is spontaneously broken for 
-3 < c 3 < -1. 
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